\ 


I 


I 


DYNAMICAL AND PHOTOMETRIC INVESTIGATION 
OF COHETARY TYPE II TAILS 

Grant NGR 09-015-159 

Semiannual Progress Report No. 6 
For the period March 15 to September 14, 1974 


Principal Investigator 
Dr. Zdenek Sekanlna 


(NA3A-CR-14C445) uYNA: 1ICAL AND 
PdOTCiITThlC INViSi'IGATION OF CO«ErAFY TYPE 
2 TAILS ScB-idniiUdl Pro»jreso Lu^poit, 15 
tiJLir* (Sit. iLvh SOI; \s^ro^hysi.C;il. 

Observatory) 36 p iiC GSCL ■'.3 A G3/3C 


Prepared for 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


N74-34294 


’ I 1 » . ' 

•« Vo ^ 

5v4i: 



Smithsonian Institution 
Astrophyslcal Observatory 
Cambridge, Massachusetts 02138 



DOMICAL AMD PHOTOMETRIC INVESTIGATION 


OF CGMETARY TYPE II TAILS 

Grant NGR 09-015-159 

Semiannual Progress Report No. 6 
For the period March 15 to September 14, 1974 


Principal Investigator 
Dr. Zdenek Sekanlna 


Pr''*^ared for 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


Smithsonian Institution 
Astrophysical Observatory 
Cambridge, Massachusetts 02138 


j 



1 


1 


TABLE OF CONTENTS 

ABSTRACT iii 

PART A. COMPARISON OF THE WORKING MODEL FOR THE ANTITAIL OF COMET 

KOHOUTEK (1973f) WITH GROUND-BASED PHOTOGRAPHIC OBSERVATIONS . 1 

!• Introduction 1 

II. The Cerro Tololo photographs 1 

III. Photographic photometry of the Cerro Tololo plates. The 

technique 4 

IV. Photographic photometry of the Cerro Tololo plates. The 

results 9 

V. Calibration stars for the Cerro Tololo plates 15 

VI. Preliminary physical interpretation of the observed radial 

profiles of the antitail 20 

VII. References 25 

PART B. OTHER ACTIVITIES IN THE REPORTED PERIOD 26 

P.yiT C. PLANS «'0R THE NEXT PERIOD (SEPTEMBER 15, 1974 TO MARCH 14, 

1975) 27 

Appendix A . Effects of differential refraction on the distance and 

directional determinations in the antitail at low altitudes . . A-1 

Appendix B . Photometric calibration of the sky background B-1 

Appendix C . Normalized size-density distribution of vaporizing par- 
ticles C-1 

ii 


1 




ABSTRACT 

Photographic observations of the antltall of Comet Kohoutek (1973f ) • 
obtained by F. D. Miller, University of Michigan, at the Cerro Tololo 
Interamerican Observatory, are photometrically studied and the results 
compared with the working model of the antltall, formulated by the writer. 
The applied technique of photometric reduction is described and the radial 
and transverse profiles of the antltall, corrected for the effects of the 
variable intensity of the sky background, are derived. The most important 
result reached so far is a quantitative confirmation of the previously 
suggested hypothesis, arguing that dust particles in the antitail suffered 
a significant loss in radius due to evaporation near the perihelion passage. 
We find that only particles initially larger than 0.1 to 0.15 mm in diameter 

-3 

(at an assumed density 1 g cm ) survived. Numerically, however, this 
result is still tentative, because the dynamical effect exerted by particle 
evaporation remains to be accounted for. 

Other activities undertaken in the reported period are also mentioned, 
and the plans for the next semiannual period are briefly listed. 
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PART A. C(»1PAR1S0N OF THE WORKING MODEL FOR THE ANTITAIL OF CCMET W)HOUTEK 


(1973£) WITH GR(XJND-BASED PHOTOGRAPHIC OBSERVATIONS 

I. Introduction 

A working model has recently been formulated for the antltall of 
Comet Kohoutek (Sekanlna 1974 > referred to hereafter as Paper I; see also 
Semiannual Report No. 5). The model Is based on the Flnson-Probsteln (1968) 
theory of dust comets and fits reasonably well the semlquantltatlve descriptions 
of the antltall by a number of observers. Yet some fundamental tests of the 
working model of the antltall, such as radial brightness variations, transverse 
profiles and time variations In the surface brightness, were missing In Paper I 
because quantitative photometry was unavailable at the time of the paper's 
completion. Some of these tests are Included In the present report. 

The photographic material studied here consists of the plates secured 
at the Cerro Tololo Interamerlcan Observatory, near La Serena, Chile, under 
the su 'ervlslon of F. D. Miller. 

II. Tne Cerro Tololo photographs 

A total of 12 plates were taken at Cerro Tololo with the University of 
Michigan Curt Is- Schmidt telescope (61/91 cm, f/3.5, 96V6 per millimeter) 
between January 16 and February 15, 1974. Five plates, taken on January 23, 

24, 26 and on February 12 and 15, are 098-02 panchromatic emulsions combined 

O 

with an RGl filter to provide the maximum sensitivity near 660QA. The other 
seven plates, taken on January 16, 17, 20, 21, 23, 25 and 27, are standard 
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103a0 emulsions without any filter. All the plates were calibrated with 
the use of 14 photometric spots from the Weston tube sensltometer. 

The red-sensitive plates have been of primary interest, both because 
of their spanning a longer Interval of time and because of their better 
showing the antltall. The guiding of the January 24 plate was somewhat 
Imperfect, so that only the remaining four plates - two in January and two 
in February - have been photometrically studied. The four observations are 
listed in Table I. 

Of much concern is of course the very low altitude of the comet. At 
the end of the exposure, the comet was only slightly over 10^ above the 
horizon on plates 15520, 15687 and 15703, and about 6?5 above the horizon 
on 15554! Only thanks to the excellent sky at Cerro Tololo the plates are 
still photometrically valuable. Yet the very low altitude of the comet is 
a source of two potential complications for the photometry of the plates, 
namely the effects of differential refraction and strongly variable sky 
brightness. 

The effect of differential refraction has been quantitatively analyzed 
following the procedure described in Appendix A. Numerically, refraction has 
been found to affect the angular distance from the nucleus on the four plates 
by less than 1%, and the positional angle by not more than 0?1 at the angular 
distance of I*’ from the nucleus and not more than 0?2 at 2°. 

The effect of variable sky brightness proved considerably more severe. 

It so happened that the brightest section of the antltall directed almost 
exactly toward the horizon on all the plates, so that the maximum brightness 
gradient of the sky projected unforeshortened along the antltall. In addition. 
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Table I. Cerro Tololo red-sensitive photographs of Comet Kohoutek 
(098-02 emulsion with an RGl filter) 
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all exposures began still during thi astronomical twilight (with the 
sun's depression between 15^ and 17**), and almost the whole 15520 exposure 
was taken during the twilight (the sun's depression at midexposure 17°). 

To demonstrate the extent of Interference, we mention that the relative 
gradient of the sky brightness on plates 15520 and 15554 reached a record 
rate of 7 per cent per degree of altitude, while at the same time the 
Intrinsic surface brightness of the antltall amounted to only some 4 > 6% 
of the sky brightness at 15 arcmln from the nucleus and about 1 - 2Z at 25 
arcmln. 

III. Photographic photometry of the Cerro Tololo plates. The technique 

The technique of reduction of the Cerro Tololo plates has essentially 
been dictated by the troublesome observing conditions. Because of the 
complicated character of the sky-brightness variations near the horizon, it 
has been virtually Impossible to use the standard two-dimensional scans to 
map the background "noise" over the area occupied by the antltall (see Fig. 1). 
Preferable - because considerably more revealing - are one-dimensional radial 
tracings, passing through the comet's nucleus and effectively covering the 
whole antltall. Each radial scan Is defined by the positional angle, and the 
"noise" variations in the respective section of the antltall can reasonably be 
Interpolated from the known background outside the antltall (see Fig. 2). 
Nearby stare of known positions have been used to fix the positional angles. 
The position of the comet on the plates was measured relative to the nearby 
stars and reduced routinely with the use of the method of dependences. The 
result has also been checked by computing the expected position of the comet 
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from Marsden's (1974) definitive orbit. The measured and computed positions 
agree within 2 arcsec on plate 15520, 7 arcsec on 15554, 8 arcsec on 15687 
and 13 arcsec on 15703. These residuals can entirely be accounted for by 
the uncertainty in the location on the plates of the nucleus in the strongly 
overexposed coma of almost 2 arcmin in diameter. The extended image of the 
comet's head is responsible for an uncertainty in the distance scale of the 
scans, amounting to several per cent at 10 arcmin from the nucleus and dropping 
down to somewhat more than 1 per cent at 30 arcmin. Together with the images 
of the orientation-fixing stars, which appear as trails a few arcmin in 
length on the four plates, the comet's image is also responsible for an error 
in the positional angle of the scans, which is estimated at not more than 
about ±0?5. Recalling the effects Introduced by differential refraction 
(Section II), we may conclude that they are much smaller than the operating 
errors and can therefore be safely neglected. 

A practically important parameter of any photometric procedure is the 
size of the scanning diaphragm. Its choice is governed by a balanced compromise 
between the requirement of a sufficiently high resolution power on the one 
hand and that of a smooth background image reasonably free from the effects 
of emulsion granulation on the other. After some experimenting, the diaphragm 
used in the radial scans was fixed at 60 by 60 microns, l.e., 5.8 by 5.8 arcsec. 

The photoelectrlcally calibrated spots have been measured with a 
Macbeth densitometer. The obtained characteristic curves are plotted in 
Fig. 3. 
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IV. Photographic photometry of the Cerro Tololo plates. The results 

The photometric procedure described In Section III has been applied to 
the four Cerro Tololo plates and the results are listed In Tables II to V. 

The smoothed brightness data are missing for areas of the antitail where Its 
photographic-density curve is very poorly determined, either because of the 
excessive faintness of the antitail or due to a strong in*" -rference with a 
star's trail. 

The data of Tables II to V have been first compared with the radial 
profiles of the working model (Fig. 4D of Paper 1). All four plates 
consistently show the effect displayed in Fig. 4 for one of the scans of 
Table II; the observed radial gradient progressively increases with distance 
from the nucleus. The surface brightness varies approximately in inverse 
proportion to the distance from the nucleus at angular distances less than 
about 10', and approaches the model's inverse 2.8 power law at about a half 
of a degree from the nucleus. The quant itiitive dependence of the radial 
gradient on the angular distance varies somewhat with the positional angle 
and to some extent also changes from plate to plate. Since the adopted 
assumption of negligibly small eniission velocities implies a close correlation 
between the radial variation of the brightness in the antltall and the 
Dopulatlon index s of the particle mass distribution m”® dm, the observed 
character of the radial gradient suggests a strongly variable s in the 
particle-size range under consideration, a case that has not been discussed 
in Paper 1. 

Although the variability of the radial gradient makes further comparison 
of the working model with our photometric results rather difficult, we have 
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Table II. Surface brightness of the antitail on plate 15520 
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used some of the data of Tables II to V to coapile Figure 5. It shows the 
profiles of the antitail at 20 arcmin from the nucleus, as a function of 
positional angle. The general impression from the comparison of the 
observations with the working model is that the antitail appears to be, 
primarily in February, broader than anticipated. When an attempt has been 
made to force the 2.8 power law through the brightness curve of the antitail 
near 30 arcmin from the nucleus, the agreement with the working model's 
transverse profiles became substantially better. However, since the bright- 
ness data at such large distances from the nucleus are rather uncertain, this 
experiment has little if any practical significance. 

V. Calibration stars for the Cerro Tololo plates 

To calibrate the surface brightness of the antltall in absolute units, 
two to three stars have been selected on each photograph. Because of very 
dense cores of star trails on Curtls-Schmldt plates, only very faint stars, 
of estimated red magnitudes 13 to 14, can be chosen for calibration. Pursuing 
a two-step procedure, we first intend to use these stars to calibrate the sky 
background (at the particular spot on the plate) and then to use the sky 
background to calibrate the antitail by converting the relative intensities 
of Tables II to V to absolute units. 

Strictly, each calibration star should have been scanned across as well 
as along its trail. However, inspection of the magnified trails of the 
calibration stars in a measuring machine has shown that virtually nothing can 
be gained by scanning the stars along their trails over simply measuring their 
lengths, once the cross scans are available. 
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RELATIVE SURFACE BRIGHTNESS (MAGNITUDES) 


1 


1 


I 


409-065 



250" 260» 270“ 280“ 290“ 300“ 310“ 320“ 


POSITIONAL ANGLE FROM NUCLEUS 

Figure 5. Transverse profiles of the antltall on the Cerro Tololo plates 
at an angular distance of 20 arcmln from the nucleus. Solid 
circles: measurements; curves: working model. 
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To get a meaningful profile of the calibration stars in the direction 
perpendicular to the trails a slit diaphragm has been used in the cross scans. 

The width of the slit in the scanning direction has been fixed at 10 microns 
(or just below 1 arcsec), a size significantly smaller than the apparent width 
of the trails of the calibration stars (which is typically 50 microns). To 
avoid the effect of emulsion granulation, the length of the diaphragm along 
the trail has been kept at 200 microns (or about 19.3 arcsec), which represents 
about 10 per cent of an average trail length on the plates. Due to differential 
refraction, the star trails are not quite straight. Each star was scanned near 
the middle of its trail and in the direction perpendicular to the tangent to 
the trail at that point. Differential refraction also affects the measured 
length of the star trails, which depend, in addition, on the star brightness. 

An example of a scan of a calibration star is shown in Fig. 6. The photometric 
reduction of the scans of this type is described in Appendix B and the obtained 
results are listed in Table VI. The coordinates of the calibration stars - 
none of which is listed in the SAG catalogue - have been determined by the 
method of dependences (using up to 10 ambient reference stars) and should be 
good to within 1 arcsec in each coordinate, except perhaps for the stars on 
plates 15554, where the accuracy may be somewhat lower. Since the magnitudes of 
the calibration stars remain to be measured, the calibration cannot be pursued 
tiii the end at present. The last column of Table VI gives the difference 

between the magnitude per square arcsec Hg of the sky background of relative 
surface brightness Ig and the magnitude of the calibration star H^. Both 
magnitudes should of course refer to the sensitivity profile of the plates with 

o 

the maximum near 66U0A. 
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VI. Preliminary physical Interpretation of the observed radial 

profiles of the antltall 

The progressively increasing gradient of the radial scans Is so far 
the best established property of the antltall. The physical significance 
of the gradient variations emerges from a relation among the modified surface 
density of the emitted particles (which Is proportional to the relative 
brightness), their emission flux and size distribution, and the area they 
occupy In the plane of the sky, as formulated by Flnson and Frobsteln (1968) 
for small particle-ejection velocities. As long as we consider only large 
particles, for which the acceleration ratio 1-p of solar radiation pressure 
to solar gravity does not exceed 'vO.Ol, the relative brightness In a particular 
direction away from the nucleus can be described by a very simple expression. 
First, the emission flux Is In this case constant, because the radial direction 
is practically Identical with a synchrone. Second, the projected area, occupied 
by the particles Is proportional to 1-w. Third, for a constant particle density 
the particle diameter 6 Is Inversely proportional to 1-y, because the 
scattering efficiency of large particles for radiation pressure Is practically 
Independent of ybe particle size. And fourth, the distance R of the large 
particles from the nucleus Is proportional to 1-y. The relative brightness 
at distance R Is then 

KR)^- f(l-y) (1-y)-^ (1) 

where f(l-y) Is related to the particle-size distribution function g(6) 
thus: 

f(l-y)'o6‘» g(6)'v-(l-y)"**g(6). (2) 
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We established in Paper I that there was a relative excesa of heavier 
particles In the antitail, which we interpreted as possible evidence for a 
significant loss of radius of the dust particles due to evaporation near the 
perihelion passage. With a sufficient body of observational information we 
now can test this hypothesis more quantitatively. 

Consider an initial (pre-evaporation) distribution function of particle 
diameters of the form 

g(6)d6 % 6"«d6, (3) 

where u is a constant. Consider further that, as a result of intense solar 
heating over a limited period of time, an outer layer %d6 in thickness 
evaporates away from each particle. This amount is indeed independent of 
the particle size for material of given structure and composition. A particle 
of initial diameter 6 has now a smaller diameter y * and the particle- 

size distribution (3) changes to 

g(Y)dy (Y+A6)““dY. (4) 

The logarithmic slope of this distribution, 

5 logo g(Y) 

t(Y) , (5) 

3 logg Y 

depends on diameter y and is related to the initial exponent u by 

t(Y) ■ 7» • (6) 

i.f 

Since Y ■ C(l-p)”^ (C is a constant and 1-u now refers to y rather than to 6), 
we can write equation (6) as 
1 1 A6 

- + (1-y). (7) 

t u Cu 
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On the ocher hand, ve now have I g(Y)« so that the radial brightness 
gradient la 

3 logg I(R) 3 log^ I(Y) 3 log. Y 

V _ . . 5 . (8) 

3 logg R 3 logg Y 3 logg R 

Deriving t(Y) for a number of angular distances on each of the scanned 
profiles, assigning the proper value of l-y (depending not only on the 
distance R, but also on the tine of observation and on the positional angle 
of the scan), and plotting the resulting pairs 1-p vs. , we can test the 
hypothesis of particle evaporation by comparing such a plot with the llneot 
relation predicted by equation (~). The plot, based on the four plates ar.u 
Including only the sections of the antltall where its relative brightness 
exceeded 2 per cent of the sky background (to avoid unrealistic gradients 
from very faint areas). Is exhibited In Figure 7. A least-square solution 
Indicates a complete absence of nonlinearity (the mean error of the quadratic 
term came out A. 5 times as large as the term Itself), and provides the 
numerical values of two Important characteristics of the population of large 
dust particles In Comet Kohoutek: exponent u of the Initial particle-size 

distribution [eq. (3)] and the evaporation loss In particle diameter A6. The 
two quantities derived from both the separate plates and the whole set are 
listed In Table VII. The table also Includes the population index s of the 
differential particle-mass distribution m'^dm, which relates to the partlcle- 
slze exponent u of (3), 
u + 2 


and the acceleration ratio corresponding to the largest particles 

that evaporated completely (i.e. , whose 6 ■ A5) . 
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Figure 7. The plot of the particle acceleration 1-u versus the logarithmic gradient 
t, of the size distribution of the vaporizing i^articles in the antitail. 
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Taken at their face value, the data of Table VII appear to Indicate 
rather consistently that the population Index Is near 2.3 and the size- 
distribution exponent u near 5, and that only the particles Initially 
larger than 100 to 150 microns In diameter survived the exposure to solar 
heat. Comparison with Figure 7 of Paper I suggests that the effective 
vaporization heat of the antltall particles amounts to about 46 kcal/mole. 
However, there Is an uncertainty Involved In these results, which Is due 
to the fact that any particle evaporation also affects the particle's future 
dynamics, because It Implies a change In the magnitude of the radiation 
pressure. Therefore the present results are still preliminary, and the 
derived loss rate of particle radius due to evaporation may only be an 
order-of-magnitude estimate. In the Improved model of the antltall the 
dynamical effect of particle evaporation (variable 1-y) must be accounted 
for. 
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PART B. OTHER ACTIVITIES IN THE REPORTED PERIOD 


Besides the extensive investigation of the antitail of Comet Kohoutek, 
the following activities were also undertaken. 

1. Computer programs 

The package of computer programs, used for the calculation of particle 
orbits, now works for particle accelerations 1-p exceeding 1 (negative 

gravitation), and for both short-period and nearly-parabollc orbits of 
comets. 

2. A general study of the anoBialous tails of comets 

This study, continued from the previous period, has resulted in a 
paper "The prediction of anomalous tails of comets", published in Sky and 
Telescope , vol. 47, pp. 374-377. 

3. The anomalous tails of short-period comets 

The general rules for the appearance of the anomalous tails have been 
applied to all the short-period comets of the past; many of these comets should 
have displayed an anomalous tail on a number of occasions, but no observations 
of such tails have ever been reported. It is tentatively concluded that the 
short-period comets fail to emit relatively heavy particles (of a submillimeter 
and larger size) in significant amounts. 

4. The problem of the tails of the distant comets 

The results of the hypothesis of the icy tails of the distant comets, 
proposed in Semiannual Progress Report No. 4 (March 15 to September 14, 1973), 
has been updated and will be shortly submitted for publication in Icarus under 
the title "A Study of the Icy Tails of the Distant Comets". 
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PART C. PLANS FOR THE NEXT PERIOD (SEPTEMBER 15, 1974 TO MARCH 14, 1975) 


1. To continue (and complete as soon as possible) the study of the anti- 
tail of Comet Kohoutek. 

2. To study the antltall of Comet Arend-Roland (1957 III) and to compare 
the results with those for Comet Kohoutek. Three high-quality photo- 
metrically calibrated plates of 1957 III are available from the 
collection of Dr. F. D. Miller, University of Michigan. 

3. To study the motions of dust particles under the effects of evaporation 
(variable 1-y) , and to derive the vaporization parameters (vaporization 
rate; latent heat of vaporization of the material). 

4. To study the synchronic bands In the tails of some comets in terms of 
the hypothesis of vaporizing particles. 

5. To complete (if time allows) the statistical study of the anomalous 
tails of comets, aimed at our better understanding of the process of 
ejection of heavy particles from comets. 

It is not believed that all the above problems will be completely 
solved by the end of the next semiannual period. The actual priority to 
each of the above points of investigation will be assigned during the 
period so that the whole project is run most efficiently. 

In addition, at the invitation of Dr. B. D. Donn, Goddard Space Flight 
Center, a review paper will be given on the progress in the study of come- 
tary dust tails at the lAU Colloquium No. 25 (Study of Comets), held in 
late October 1974. 
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Appendix A « Effects of differential refraction on the distance and 
directional determinations in the antltall at low altitudes 

Because of the sharp increase in the refraction correction of the 
spherical coordinates at low altitudes. It Is desirable to establish the 
effects of refraction on the distance and directional scales In the antitail. 
Let 5^ be the refraction-free coordinates of the comet's nucleus, 

I t t I 

® 2 * ^2 ® point in the antltall, anci r.i. tlarly a^, 6 ^, a^, 6 ^ 

corresponding coordinates affected by refraction. The corrections fot 
refractior are therefore Aa^ ■ “ ^i~^i “ 1»2). Since we 

consider only small arcs, we can write for the refraction-free angular distance 
A of the point (o^, 62 ^ from the nucleus 

A^ * [(a^+Aa^)-(a^+Aa^)]2cos^*s[(Aj+A5^)+(6^+A62)]+[(A*+A62)-(6*+A6^)]^. (A-1) 

Denoting the refraction-affected distance of point (o^, 62 ^ ®® 
realizing that 

cos ^[( 5 j+A 6 ^)+( 6 ^+A 62 )]** cos cos sj, 

cos^ b(A|+ 6 *)+( 62 “^X^^ * 

|Aoi^-Aa |«|a^_a|l, (A-2) 

1A62-A6J«16;-6;1, 

A2_^'2fts 2A' (A-A') 

we find the refraction correction in distance to amount to 

A-A' i.^[(o*-a|)(Aa2-Ao^)cos26|+(6^-fi|)(A6^-A6^)]. (A-3) 

Similarly, the refraction-free positional angle p of the point ( 02 ,^ 2 ) relative 
;o the nucleus is 

tan p ■ (c^+Aa^) Js[(A^+A 6 ^)+( 62 +A 62 ) 3 . (A-4) 

(62+A62)"(52+A6j) 


A-1 



Because of (A-2), its relation to the refraction-affected positional 
angle p' is expressed as 


tan p ■ tan p' 


(1 + 


Aa^-Aoi 



A62-A6 



(A-5) 


and, with p and p' in degrees, the refraction correction in direction 
becomes 


p-p' - 28?65 sin 2p' ( 


Aa -Aa 
2__L 

I I 
®2"“l 


AA^-AAj 

6*-6' 

2 1 


(A-6) 


The corrections Aa^, A5j^ follow from the classical equations of the 
spherical astronomy: 

Aoj^ = -cos sin t^ sec 6^ cosec 2 ^ 

A6j = (tan 6i cot - sin <|) sec 6^ cosec z|)R^, (A-7) 


where (|i is the geographic latitude of the observing place, t^ and z^ the 
hour angle and zenith distance, both refraction-affected, and R the refraction 
correction in the zenith distance, which, except very near the horizon, can 
be approximated by 

R = A tan z' + B tan^ z', (A-8) 

where A and B are constants. 
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Appendix B . Photometric calibration of the sky background 


We assume that a star's Image on a calibrated plate, guided on a moving 
object (comet), shows up In the form of a trail of length L and characteristic 
width W (arcsec). We further assume that a photometric scan across a section 
of the trail Is available and that the diaphragm used has been a silt of 
length t and width w (arcsec), where t«L and w<<W. Converting the photo- 
graphic densities across the measured profile of the star to relative Intensities 
and fixing the level of the sky-background noise on the scan, we wish to 
estimate the contribution of the scanned section of the star to the brightness 
of the ambient background Ig. We assume that the excess relative Intensity 


I across the star's profile fits a probability-type curve, 

I(x) = I* exp [- ], 

2o2 


(B-1) 


where x (arcsec) is read from an arbitrary origin x , I is the peak relative 

o it 

brightness at the center of the star's profile, and a Is the width factor 
(arcsec). Establishing I(x) for a number of x from the scan we can determine 
and a by least squares. Equation (B-1) can be written as 

(B-2) 


loggl(x) » A + Bx + Cx^, 


where 


loggl* - 


2a‘ 


C - - 


so that 


20-^ 


exp (A - 


4C 


(B-3) 


(B-4) 


B-1 





.1 


] 


1 


I 


and 

a - (-20“**. (B-5) 

The Integrated intensity in the scanned section of the star's trail is 


and 

the 


I dx - (2ir)**I^o, (B-6) 

the brightness contribution of the whole trail of the star Is related to 

star's magnitude H : 

* 


(2ir)**I o -j- - c. lO"®-"***, 


(B-7) 


where c is a constant. The surface brightness of the ambient sky background 
Hg (magnitude per square arcsec) can similarly be expressed through the 
relative brightness and the slit dimensions 

— ^5— = c. 10“®*"®B. (B-8) 

tv 


The elimination of c from (B-7) and (B-8) then gives immediately the result: 


Hg - + 2.5 log^jj [(2ir)*5 I^aLw/lg]. 


(B-9) 
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Appendix C . Nonaallzed size-density distribution of vaporizing particles 

Consider particles, whose diameters 6 and densities p have originally 
satisfied a distribution law 

g(p6) d(p6) - Cg(p6)““ d(p6), (C-1) 

where u and Cg are constants. Consider further that the particles vaporize 
until their final diameters y and densities o are related to the pre-evapo- 
ration values thus: 

oy • p6 - A(p6), (C-2) 

where A(p5) is a constant. Obviously, all the particles for which p5 < A(p5) 
vaporize completely, whereas the post-evaporation distribution of the 


particles with pA > A(p6) becomes 

g(OY) d(oy) » Cg [oY + A(p6)]““ d(cTY) . (C-3) 

In order that this distribution be normalized, 

7 g(OY) d(0Y) = 1, (C-4) 

o 

the constant Cg must satisfy the relation 

Cg » (u - 1) [A(p6)]“-1. (C-5) 

The distribution of radiation-pressure accelerations 1-p, associated 
with the particle size-density distribution (C-3) is 

f(l-u) d(l-w) ■ C^ (oy)** g(OY) d(l-p), (C-6) 

or because ay « C(l-p)“^ 

f(l-M) d(l-w) - Cf . [1 + d(l-M), (C-7) 


where C^ “ ^f*^g (l-u) mav radiation-pressure acceleration exerted 

on the largest particles that later evaporated completely [ p6 = A(p6)]. 
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1 


1 


] 
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Normalizing distribution (C**7) we obtain 

Of - 3 (l-u)^ • C“-*, (C-8) 

SO that 


f(l-y) d(l-y) - 3 ('*3^) (l-y)3^ [1 +-(xES) (C-9) 

Finally, the important parameter appearing In the expression for the 
mass-loss rate of the dust and numerically equal to the ratio of the mean 
cube to the mean square of the size-density distribution [see eqs. (8) and 
(10) of Paper l], becomes 




3C 

(u-4) (l-y)niax 


3 A(p6) 
u-A 


(C-10) 
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